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Superoxide dismutase (SOD) is a critical enzyme associated with controlling oxygen toxicity arising out of 
oxidative stress in any living system. A hyper-thermostable SOD isolated from a polyextremophile higher 
plant Potentilla atrosanguinea Lodd. var. argyrophylla (Wall, ex Lehm.) was engineered by mutation of a 
single amino acid that enhanced the thermostability of the enzyme to twofold. The engineered enzyme was 
functional from sub-zero temperature to >50°C, tolerated autoclaving (heating at 121 °C, at a pressure of 
1.1 kg per square cm for 20 min) and was resistant to proteolysis. The present work is the first example to 
enhance the thermostability of a hyper-thermostable protein and has potential to application to other 
proteins for enhancing thermostability. 



Copper, zinc superoxide dismutase (Cu,Zn SOD; EC 1.15.1.1) catalyzes a two-step dismutation of the toxic 
superoxide radical (O 2 0 to molecular oxygen and hydrogen peroxide through alternate reduction and 
oxidation of copper ion at diffusion-limited catalytic rate 13 . Cu,Zn SOD forms the first line of defense in 
all aerobic organisms against reactive oxygen radicals and is vital to the survival of cells 4 . Therefore, SOD is an 
important enzyme for aerobic organisms including plants and humans and also in those applications which 
require scavenging of 0 2 5 ~ 7 . 

Eukaryotic Cu,Zn SODs display conserved three dimensional fold built on a flattened Greek key motif and are 
homodimeric with two subunits held together by hydrophobic interactions 8,9 . Several reports suggested dimeric 
structure of Cu,Zn SOD to be more stable with superior activity 1012 , though in Escherichia coli this did not appear 
to be the case 13 . Cu,Zn SODs are known stable enzymes since these tolerate 8 M urea 14 and 4% sodium dodecyl 
sulfate 15 . Various factors responsible for imparting stability/thermo-stability include relatively small solvent- 
exposed surface area, increased packing density that reduces cavities in the hydrophobic core, an increase in core 
hydrophobicity, hydrogen bonds between polar residues, increased occurrence of proline residues in loops, 
decreased occurrence of fhermolabile residues, increased helical content, increased polar surface area, extensive 
ion-pair networks and better salt bridge formation 8,1617 . 

A unique Cu,Zn SOD (WT) in P. atrosanguinea was identified by our group 5 that was functional at sub-zero 
temperature to >50°C and also the protein retained activity after autoclaving (heating at 121°C, at a pressure of 
1.1 kg per square cm for 20 min). The enzyme was identified as a dimer consisting of two asymmetric subunits A 
and B 8 . Each subunit of the enzyme consisted of eight stranded |3- sandwich (la-8h) connected by seven loops (T 
VII), one short helix and one 3 10 helix. Since increased thermostability is an important desirable feature for Cu,Zn 
SOD 6 , there is a need to develop methods to improve thermostability of the enzyme. 

Several studies showed the importance of single amino acid substitution in modulating kinetic properties of 
enzymes 11,18-20 . We developed seven mutants by replacing amino acids at targeted positions by site-directed 
mutagenesis wherein both subunits of the enzyme had desired mutations. For example, amino acids in the 
electrostatic loop are important as these determine the shape and strength of the electrostatic field around the 
active site 8 . Accordingly, two important charged residues in the electrostatic loop, Leu-132 and Ser-135 were 
replaced by Glu and Lys in two separate mutations. Alanine at position 4 is yet an important amino acid, wherein 
mutation is reported to cause amyotrophic lateral sclerosis in human 21 . In WT, glycine is present at position 4 and 
this was mutated to Ala and He in two separate mutations. All of these mutations decreased the specific activity of 
the enzyme and also the tolerance to autoclavability got reduced (Supplementary Fig. SI). 
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Presence of cysteine in SOD has drawn considerable attention. 
Mutation of cysteines involved in disulphide linkage results in desta- 
bilization of Cu,Zn SOD 22 . However, mutation of free cysteine 
increased greater resistance to thermal inactivation 18 . Experimental 
evidence for the enhanced thermo-stability of the Cys (free) — > Ala 
substitution of enzymes has been reported for Cu,Zn SOD from Bos 
taurus 1 *, Danio rario 19 , Xenopus laevis 20 , and an esterase from 
Bacillus stearothermophilus 2 * . Cu,Zn SODs from yeast and E. coli, 
which have no free cysteine residues, show partial and full revers- 
ibility after heating at temperatures above their melting tempera- 
tures 24,25 . However, Cu,Zn SOD from Photobacterium leiognathi is 
denatured irreversibly despite having no free cysteine residues 26 . 
These results suggested that the free cysteine might or might not 
be related to the thermal stability of Cu,Zn SOD. Further, studies 
for characterization of thermo-stability of Cu,Zn SODs had been 
carried out with emphasis on thermodynamic parameters and the 
absence of data on enzyme activity posed a question on functionality 
of the enzyme. 

WT enzyme had three cysteine residues at positions 56, 95 and 
145, out of which Cys-56 and Cys-145 participated in the formation 
of disulfide linkage, whereas Cys-95 remained free 8 . The role of free 
cysteine in regulating the kinetic properties varies, possibly depend- 
ing upon structure of the protein and hence the amino acid sequence 
of the protein. So far structural information of Cu,Zn SODs from 
three plant species namely P. atrosanguinea (PDB ID: 2Q2L), 
Spinacea oleracea (PDB ID: 1SRD), and Solarium lycopersicum 
(PDB IDs: 3 MKG, 3 KM1, 3 KM2, 3 HOG) is known. However, 
the study of molecular determinants for structural-functional aspects 
of plant SODs is not yet reported. 

The objective of the present work was to study the role of cysteine 
in governing the kinetic properties of Cu,Zn SOD of P. atrosanguinea 
with the background that there is no such knowledge existing for 
plant Cu,Zn SODs, and that plant and the animal Cu,Zn SODs are 
structurally very different to each other 7 . Since SODs have enormous 
implications to the living world due to their role in scavenging super- 
oxide radicals, the present work provides a highly flexible yet stable 
Cu,Zn SOD at varying temperature. 

Results 

Native SOD enzyme is being referred to as wild- type (WT), whereas 
mutant SODs in which Cys-56, Cys-95, Cys-145 were replaced with 
Ala will be referred to as C56A, C95A, and C145A, respectively in the 
present work. 

The free cysteine on the interface potentially mediates key dimer 
interactions. To understand the possible role of free Cys-95, the 
interface residues and their interactions were studied in silico. Using 
relative difference in solvent accessibility criterion of POLYVIEW 
server 27 , residues on the dimer interface were identified. Eight 
residues on chain A (Thr-17, Leu-19, Thr-31, Asn-33, Asp-91, Thr- 
93, Cys-95, Phe-97) and 6 residues (Glu-12, Asn-54, Met-57, Leu- 132, 
Gly-140, Gly-141) on chain B were predicted to be on the dimer 
interface (Supplementary Fig. S2). It was observed that the dimer 
interface has 8 hydrophilic residues (out of 14 predicted to be on 
the interface) including free Cys-95. Graph-theoretical analysis 
indicated that Cys-95 residue on the interface is in close structural 
proximity of other interface residues via short- and long-range 
contacts. Cys-95 was in direct long-range contact with two other 
residues on the interface: Thr-31 and Asn-33. Cys-95 interacted 
via back-bone contacts with Thr-93, and through non-covalent 
interaction with Phe-97 (one-residue-apart) and Asn-91 (three- 
residue-apart). These contact network based observations, relying on 
structurally critical long-range contacts, proposed a role of Cys-95 in 
oligomerization of WT through interface forming and stabilizing 
contacts of WT dimer (Supplementary Fig. S3). 



The cysteines involved in native disulphide bond formation are 
part of 'Critical Interaction Residue Clique (CIRC)'. A CIRC was 
identified by graph-theoretical modeling of the WT native-state 
structure. The long-range interaction network (LIN) model suggested 
a CIRC comprising of four residues that formed a strong long-range 
clique. Residues Cys-56, Val-117, Ala- 144, and Cys-145 were predicted 
to be critical for the structural stability (Supplementary Fig. S4). Thus, 
graph-theoretical modeling located the cluster of residues that holds 
different secondary structures together, apart from identifying the 
native disulphide forming cysteines (Cys-56 and Cys-145) as a part 
of CIRC. 

Mutation of free cysteine enhances monomer-dimer ratio. WT 

and C95A existed in monomeric and dimeric form when electro- 
phoresed on 15% SDS-PAGE under non- reducing conditions 
[without 2-mercaptoethanol (2-ME)]. Whereas, C56A and C145A 
exhibited multiple bands due to possible oligomerization/aggrega- 
tion. A noticeable result was a substantial increase in monomer to 
dimer ratio for C95A as compared to the WT (Fig. la). Monomer and 
dimer of WT and C95A exhibited SOD activity on SDS-PAGE, 
whereas no SOD activity was detected on the lanes loaded with 
C56A and C145A (Supplementary Fig. S5). Isolation of monomer 
and dimer by size exclusion chromatography (Fig. lb) followed by 
protein quantification showed that monomer to dimer ratio for the 
WT and C95A was 1.3 ± 0.2 and 42.8 ± 1.0, respectively (Fig. lb). 
Monomer and dimer of the corresponding peaks were also 
confirmed by SDS-PAGE analysis under reducing and non- 
reducing conditions (Fig. lc), and the identified peaks exhibited 
SOD activity (Fig. lb). Since C56A and C145A did not exhibit 
SOD activity, monomers, dimer and oligomers were not purified 
for these two cases. It was observed that the monomeric form of 
SOD migrated as a closely spaced doublet on SDS-PAGE. The 
closely spaced protein bands in doublets were confirmed to be 
identical as evidenced by N-terminal amino acid sequencing (des- 
cribed in Supplementary Methods) that exhibited identical amino 
acid sequences MRGSHHHHHHGS. Data was further supported by 
mass analysis using a MALDI-TOF (AXIMA-CFR Plus spectro- 
meter, Shimadzu-Kratos, UK) that yielded a mass of 17.39 kDa for 
WT (Supplementary Fig. S6). The value was in agreement with those 
obtained by SDS-PAGE (Supplementary Fig. S7) and in silico 
analysis (www.expasy.org) that deduced a mass of 17.4 kDa for the 
WT. WT/C56A/C95A/C145A is induced in a relatively higher 
quantity as compared to the individual native proteins in E. coli 
(Supplementary Fig. S7) and it appears that a part of the induced 
protein acquired a differential tertiary structure leading to dif- 
ferential levels of SDS-loading and hence different migration of 
polypeptide-SDS complex on SDS-PAGE 28 . 

Cysteine replacement affects SOD activity at different tempera- 
tures without affecting pH optima. Enzymatic activity of SOD was 
checked on a non-denaturating PAGE by activity staining followed 
by spectrophotometric-based assay. SOD activity for WT and C95A 
was detectable as opposed to C56A and C145A, which did not exhibit 
any activity (Fig. 2a). Spectrophotometric-based data showed a 
higher specific activity of C95A as compared to WT by 5.0 ± 
0.6%. After autoclaving, WT lost 17 ± 4% of the un-autoclaved 
activity, whereas this loss in activity for C95A was only 8 ± 2% 
(Fig. 2a). Since no SOD activity was detectable for C56A and 
C145A, these were not evaluated for other enzymatic studies. WT 
exhibited maximum activity at 0-4°C whereas it was at 10°C for 
C95A (Fig. 2b). WT and C95A were active at zero and sub-zero 
temperatures. Autoclaved WT exhibited 50% of maximum activity 
of un-autoclaved enzyme between 30-40°C, whereas it was observed 
between 50-65°C for C95A. 

Interestingly, assay of WT and C95A at pH ranging between 6 to 9 
showed that the mutation of Cys-95 did not affect the pH optima of 
7.8 (Supplementary Fig. S8). 
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Figure 1 | Analyses of aggregation pattern (including monomer to dimer) of WT and mutant enzymes. One jig of each WT and mutant proteins was 
loaded separately onto 15% SDS-PAGE in the presence ( + ) and absence (— ) of 2-ME in the sample loading buffer 50 (a). The name of sample is given 
above top of each panel. M, molecular weight marker procured from Fermentas, USA; the protein markers were- fi-galactosidase (11 6000 Da), bovine 
serum albumin (66 200 Da), Ovalbumin, 45 000 Da), lactate dehydrogenase (35 000 Da), REase Bsp981 (21 000), p- lactoglobulin (18 400) and lysozyme 
(14 400 Da)]. SODs analyzed as in (a) were also stained for SOD activity (Supplementary Fig. S5) wherein the said activity was detected for WT and C95A, 
but not for C56A and C145A. Monomer and dimer observed in panel "a" were separated by size exclusion chromatography (b) for which the purified 
protein after Ni-NTA chromatography was used as described in Methods section. Numerical values below monomer and dimer represent the amount of 
protein (ug) as measured by Bradford's reagent" 8 (Sigma, USA), whereas those in parentheses are the specific activity of SOD (units/mg protein). 
Separated monomer and dimer of the corresponding peaks in panel "b" were also checked on a 15% SDS-PAGE in the presence ( + ) and absence ( — ) of 
2-ME included in the sample loading buffer (c). Due to large variation in the protein content in the peaks, equal quantity of protein could not be used for 
the loading. M, molecular weight marker as in "a"; WT, wild type; C56A, cysteine at position 56 substituted with alanine; C95A, cysteine at position 95 
substituted with alanine; C145A, cysteine at position 145 substituted with alanine. 
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Figure 2 | Effect of cysteine replacement on the activity of SOD. Activity staining gel showing the effect of Cys^Ala replacement on SOD activity (a). 
Equal quantity (300 ng) of each un-autoclaved and autoclaved WT and mutants was loaded onto 10% native-PAGE and stained for activity 46,47 . Names of 
the protein samples are shown at the top of each panel whereas specific activities are given at the bottom of the panel. Enzyme activity was expressed as 
unit/mg of protein. Values were mean ± SE of three separate replicates. ND, not detectable. Effect of assay temperature on SOD activity is given in panel 
(b). Per cent relative activity with respect to the maximum activity across the temperature range is shown. Solid and dashed lines represent un-autoclaved 
and autoclaved enzymes, respectively. Values were mean ± SE of three separate replicates. WT, wild type; C56A, cysteine at position 56 substituted with 
alanine; C95A, cysteine at position 95 substituted with alanine; C145A, cysteine at position 145 substituted with alanine. 
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Mutation of free cysteine improves resistance to thermal inactiva- 
tion. Residual SOD activity of WT was 48.0 ± 1.4%, as compared to 
69.0 ± 3.4% for C95A after heating for 160 min at 80°C. Thermal 
inactivation data exhibited the most linear relationship when plotted 
at first-order kinetics (Fig. 3) which yielded correlation coefficient (r) 
of -0.97 and -0.98 for WT and C95A, respectively. Thermal 
inactivation rate constant (k^) decreased from a value of 4.6 ± 
0.2 X 10~ 3 mfrr 1 for WT to 2.4 ± 0.3 X 10~ 3 mfrr 1 for C95A. 
Concomitantly, the half-time of thermal inactivation (t 1/2 ) at 80°C 
increased from 151 ±6 min for WT to 345 ± 36 min for C95A 
(Fig. 3 and Supplementary Table SI). The kinetic data showed 
improved resistance of C95A towards thermal inactivation, as com- 
pared to the WT. 

WT and C95A tolerate limited proteolysis. WT and C95A 
exhibited proteolytic tolerance. The digestion patterns as detected 
by SDS-PAGE revealed partial hydrolysis of WT and C95A with 
trypsin, chymotrypsin and papain (Fig. 4a). Of particular interest 
was the maintenance of catalytic function by WT and C95A 
(Fig. 4b,c). WT and C95A retained >85% of activity after 
proteolysis. C56A and C145A, which did not exhibit SOD activity, 
were susceptible to proteolytic cleavage as evidenced by cleaved 
bands of low molecular weight on SDS-PAGE (Fig. 4a). 

Circular dichroic spectroscopy (CD) analysis reveals enhanced 
conformation stability in C95A. CD-spectra for WT and C95A 
had a broad minimum at 208 nm (Fig. 5a) suggesting high per- 
centage of P-sheet content as observed in human and most other 
eukaryotic SODs. For C56A and C145A, shift of CD-signal towards 
200 nm was observed eliciting enhanced randomness in the struc- 
tures. CD-spectra revealed highest conformational stability of C95A 
towards autoclaving. After autoclaving, C95A retained maximum 
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Figure 3 | Thermal inactivation kinetics of SOD. 1 .0-ml solution of 
0.2 mg/ml of each enzyme in 50 mM potassium phosphate buffer, pH 7.8, 
was incubated at 80°C for 160 min. Samples (100 ul) were removed after 
every 20 min and stored on ice. After 160 min, all samples were assayed for 
SOD activity under standard assay conditions using NBT-based assay. The 
data were fit as zero, first, and second order reactions 18 , wherein the first- 
order reaction gave the most linear relationship and shown in this Fig. 3. 
SOD activity was not detected for C56A and CI 45 A, and hence 
corresponding activity curves are not shown. Solid and dashed lines 
represent WT and C95A, respectively. Values were mean ± SE of three 
separate replicates. WT, wild type; C95A, cysteine at position 95 substituted 
with alanine. The above figure was also used to calculate various kinetic 
parameters as shown in Supplementary Table SI. The rate constant 
(min -1 ) and f 1/2 for first-order thermal inactivation were determined 51,52 
from "equation (1)" and "equation (2)", respectively (equations described 
in "Methods" section). A t , residual activity that remains after heating the 
enzyme for time t; A), initial enzyme activity before heating. 



|3-sheet content followed by WT (Supplementary Table S2). C56A 
and C145A had very low CD-signal after autoclaving and were the 
most labile structures (Fig. 5a and Supplementary Table S2). 

Thermal-stability of un-autoclaved and autoclaved WT and 
mutants was further studied by recording CD-spectra at different 
temperatures upto 100°C (Fig. 5b). WT and C95A retained the sec- 
ondary structure upto 60°C and a loss in their structure was recorded 
at higher temperatures. C56A and C145A exhibited randomness in 
structure at temperatures above 50°C. 

Cysteine replacements affect the three dimensional structure 
related parameters of SOD: an in silico modeling analysis. Based 
on the three-dimensional structure of WT, obtained from X-ray 
crystallography with 2.1 A resolution 8 , models of WT, C56A, C95A, 
C145A were constructed. The deviations between these models were 
generally small for the backbone atoms (rmsd from 0.2 to 0.6 A). The 
constructed models also showed high verification scores when 
estimated using "Verify 3D Profile" tool of the DS. Differences were 
obtained for these three mutants for DOPE score, PDF total energy, 
and total molecular volume (Supplementary Table S3). The regions of 
deviations were visible after superimposing the generated models 
using WT (2q21, subunit a). For C95A, the DOPE score and PDF 
total energy was found to be lowest amongst the four enzymes 
evaluated suggesting a more stable structure. Also, these mutants 
were analyzed using Build Mutant tool of the DS. The energy 
calculations revealed that C56A and C145A showed increase in the 
PDF total energy (33%) with respect to the WT, whereas C95A 
resulted in slight decrease (1%) in the PDF total energy. 

Discussion 

Role of cysteine in modulating kinetic properties of Cu,Zn SOD is 
still an enigma 13 ' 18-2 "' 23,24 . Enormous work has been carried out in 
animal system, but there are no reports on the role of cysteine in 
any of the Cu,Zn SODs derived from plants. Since amino acid 
sequences of Cu,Zn SODs in plant and animals vary significantly 
hence any conclusions drawn for SODs derived from animals may 
not be valid for those derived from plants. Cu,Zn SOD from P. 
atrosanguinea provides a unique opportunity to understand the role 
of free Cys-95 since the native enzyme retains more than 80% of its 
activity upon autoclaving and functions at temperature ranging 
between sub-zero temperatures to more than 50°C with maximum 
activity at 0-4°C 5,8 . WT has three cysteines, out of which Cys-56 and 
Cys-145 are involved in disulphide linkage, whereas Cys-95 remains 
free. 

The analysis on relative difference in solvent accessibility criterion 
suggested free Cys-95 to be present on the dimer interface 
(Supplementary Fig. S2). The presence of Cys-95 on dimeric inter- 
face of WT enzyme subunits is likely to play an important role in 
subunit association. Contact network based observations, based on 
structurally critical long-range contacts, suggested a possible role of 
Cys-95 in oligomerization of WT enzyme (Supplementary Fig. S3). 
Possibility also exists that Cys-95 would interact with another Cys-95 
through disulphide linkage to form oligomer. Interestingly, Cys-95 
adopts a double conformation in subunit B 8 . Increased number of 
conformation possibilities decreases the protein stability by increas- 
ing the entropy of the protein in unfolded state 29 . Thus restricting 
conformational possibilities of the unfolded state is important in 
lowering its entropy and prevent unfolding. 

The in silico based propositions on the role of Cys-95 on oligo- 
merization were evaluated by replacing free Cys-95 with Ala. Alanine 
was a choice in this replacement, since the amino acid eliminates side 
chain beyond the P carbon. Also, alanine neither imposes extreme 
electrostatic/ steric effects, nor alters the main-chain conformation 30 . 
This replacement increased monomer to dimer ratio to 33 fold in 
C95A as compared to that in WT (Fig. la,b). Interestingly, the spe- 
cific activity of C95A was higher as compared to the WT (Fig. 2a), 
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Figure 4 | Effect of incubation of SOD with proteolytic enzymes. SOD was incubated in the absence (-) and presence ( + ) of 1/20 w/w of trypsin, 
chymotrypsin or papain in 50 mM potassium phosphate buffer, pH 7.8, and 37°C for 1, 2, and 3 h and run on a SDS-PAGE (a). M, molecular weight 
marker procured from Fermentas, USA; the protein markers were- fi-galactosidase (11 6000 Da), bovine serum albumin (66 200 Da), Ovalbumin, 45 
000 Da), lactate dehydrogenase (35 000 Da), REase Bsp981 (21 000), |3- lactoglobulin (18 400) andlysozyme (14 400 Da)]. Full-length gels are shown in 
Supplementary Fig. S9a. Activity staining of protease treated SODs is shown in panel "b". Protease treated SODs were run on native rather than SDS 
containing poly-acrylamide gel. Therefore, the migration pattern of SODs shows some variation between the panels (a) and (b). SOD activity was not 
detected for C56A and C145A, and hence corresponding SOD activity stained gels are not shown. Full-length gels are shown in Supplementary Fig. S9b. 
Panel (c) depicts the effect of proteolytic enzymes on SOD activity; values are mean ± SE of three separate replicates. WT, wild type; C56A, cysteine at 
position 56 substituted with alanine; C95A, cysteine at position 95 substituted with alanine; C145A, cysteine at position 145 substituted with alanine. 



suggesting that higher monomerization favored increased catalytic 
activity. Monomer as well dimer form of the Cu,Zn SOD is known to 
retain the catalytic activity 1131 . Depending upon the species under 
consideration the catalytic activity of the two forms could vary. For 
example, dimer exhibited comparatively higher activity in Ipomoea 
batatas 31 , and Carica papaya 11 . Higher activity of dimer form was 
possibly due to improved enzyme conformation and stability of the 
enzyme 11 . Monomerization led to an abrupt decrease in the catalytic 
activity of the enzyme either by alterations in the tertiary structure 
due to extensive solvation and/or distortion of the newly exposed 
dimer interface in Cu,Zn SOD from human 1332,33 . On the contrary, 
monomeric Cu,Zn SOD of E. coli retained the full catalytic activity, 
which suggested that dimeric form is not a requirement for the 
efficient catalytic activity by Cu,Zn SOD 13 . Substitution of a charged 
amino acid in place of a hydrophobic amino acid has been found to 
be responsible for monomerization 34 . However, in C95A substi- 
tution of a uncharged polar amino acid with a hydrophobic amino 
acid (Ala) led to higher monomerization (Fig. la,b). Results were 
in consonance with the computational data that suggested the role 
of Cys-95 in interfacing with other amino acids to stabilize dimer 
of WT (Supplementary Fig. S3). The present data suggested the 
importance of higher monomer to dimer ratio for enhanced ther- 
mal tolerance as evident from the reduced (about half) k d and 
increased t 1/2 (about twofold) values for C95A as compared to 
WT (Supplementary Table SI). 

Lesser loss in activity after autoclaving and lower thermal inac- 
tivation in C95A as compared to WT suggested reduced irreversible 
denaturation in the former. The results were supported by CD data 
which showed a loss of 20% of |3-sheet content upon autoclaving WT 
as compared to the un-autoclaved control. However, C95A exhibited 
an increase of 2.5% of [3-sheet content upon autoclaving as compared 



to the un-autoclaved control (Fig. 5a and Supplementary Table S2). 
These results are suggestive of retention of secondary structure of the 
autoclaved C95A as compared to the autoclaved WT. 

Lower values of PDF physical energy and the DOPE score 
(Supplementary Table S3) also indicated C95A structure to be the 
least restrained. Since free cysteine might be the preferred site for 
oxidation of SOD molecules, it could impede the enzyme activity. It 
has been shown that human Cu,Zn SOD undergoes irreversible 
aggregation upon exposure to temperatures higher than their 
respective melting temperatures, which is attributed to the presence 
of free cysteine residues. Formation of improper disulfide-bonds, 
concomitant with cysteine oxidation has been attributed to heat- 
denatured aggregation 34 . Mutation of free Cys-95 to Ala might elim- 
inate any such possibility, and would render the protein more stable. 

It was also observed that cysteine residues of intra-disulfide bond 
(Cys-56 and Cys-145) are important for functioning of WT, since 
mutation of these residues resulted in the loss of enzymatic activity. 
The results were strengthened by CD analysis where structural loss 
was observed for these mutants (Fig. 5a,b). The loss maybe attributed 
to permanent unfolding which is usually combined with the forma- 
tion of aggregates due to the exposed hydrophobic core 35 . The 
aggregation of C56A and C145A observed on SDS-PAGE suggested 
inter or intra- molecular non-native disulfide linkages in the proteins. 
The loss in structure at higher temperatures could be assigned to the 
increased entropy with increased number of conformational pos- 
sibilities, which is usually high in the unfolded state and low in the 
folded state. It has been proposed that cysteine disulfide bridges can 
decrease the entropy of the unfolded state of a protein and as such 
acts as a stabilizing factor 36 . Disulfide bond links loop IV and (3- 
strand 8h of WT and this linkage of the secondary- structure elements 
contribute to the stabilization of the SOD fold 8 . The active site is 
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present at bottom of channel formed by two loops (IV and VII) 
which extends from )3-barrel. Mutating cysteine residue in loop IV 
might destabilize WT and thereby decreased activity was recorded by 
affecting the active site. Graph-theoretical modelling studies revealed 
Cys-56 and Cys-145 to be critical for structural stability which was 
supported by our experimental work. 

Resistance to proteolysis is a selection tool to assess protein 
stability 3738 . Also, proteolysis is used to probe energy states in proteins. 
Proteins are protected against proteolysis in their folded conforma- 
tion, unless these have unstructured regions. In order to be digested, 
the substrates must access a high-energy state where cleavage sites are 
exposed to the solvent 39 . Also proteolysis depends upon the accessible 
surface area. WT and C95A were found to be the least susceptible to 
cleavage by proteolytic enzymes trypsin and chymotrypsin at normal 
temperature. Papain partially digested both the enzymes (Fig. 4a), still 
the activity was retained upto 85 to 90% of the full activity (Fig. 4b,c). 
Results strengthened the hypothesis of a general correlation between 
the thermo- stability of proteins and their resistance to proteolysis 37 . 



Resistance of WT and C95A to proteolytic cleavage suggested that 
the protein has limited access to partially and globally unfolded 
conformations under native conditions 40 . Proteolytic resistance 
might have been imparted due to minimized occurrence of accessible 
conformations susceptible to proteolytic attack. C56A and CI 45 A 
were highly susceptible to proteolytic cleavage (Fig. 4a), possibly due 
to loss of structure of these mutants as evidenced from CD data 
(Fig. 5a,b) leading to increased access of the protein by trypsin, 
chymotrypsin and papain. 

To conclude, the present study suggested the importance of 
increased monomer to dimer ratio to enhance thermostability of 
Cu,Zn SOD, wherein mutation of a free cysteine (Cys-95) played a 
central role. Data was supported by in silico as well as other 
biochemical parameters. The engineered SOD can be used for 
developing transgenic plants tolerant to abiotic stresses 41,42 , par- 
ticularly for high temperature and drought stress, where the tem- 
perature can rise to as high as 50-60°C 43 . Herein, ROS mediated 
damages are inevitable. C95A reported in the present work, which 
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is tolerant as well as functional to higher temperature, would be a 
boon for such application. 

Methods 

Site directed mutagenesis. Coding sequence of WT (456 bp; GenBank Acc No., 
EU532614) was amplified and cloned in pQE-30 U/A expression vector 8 ' 44 . 
Polymerase chain reaction (PCR) was carried out to create mutation in SOD gene by 
using mutagenic primers (Supplementary Table S4) as described previously 45 . In 
brief, PCR mix (25 ul) contained 2.5 ul of lOx long-amplicon (LA)-Taq polymerase 
buffer, 100 nM each of dATP, dTTP, dGTP, and dCTP; 100 ng of each forward and 
reverse primers; 25 ng plasmid pQE30-UA carrying SOD gene; 0.5 ul (1.25 U) LA- 
Taq (Takara, Clontech, USA) and water to a final volume of 25 ul. PCR product was 
digested with Dpn \ and the resulting single stranded mutated DNA was transformed 
into E. coli XLl-Blue supercompetent cells (Stratagene, USA). Colonies were screened 
on LB ampicillin plates, plasmids were isolated using GenElute Plasmid Miniprep Kit 
(Sigma, USA) and sequencing was performed on an automated DNA sequencer (ABI 
Prism 3130x1, Applied Biosystems, USA) to confirm the incorporation of desired 
mutations. The mutated plasmids were transformed into E. coli M15 cells. 
Recombinant proteins were expressed and purified under native conditions by 
affinity based chromatography using Ni-NTA 44 . Since, SOD is a homodimeric 
enzyme both subunits had the desired mutations. 

SOD assay and protein estimation. Enzyme assay was performed as described 
previously 46 ' 47 except that the reaction was carried out using 1 jig of the purified 
protein per assay in a 96 well plate maintained at 4 C (optimal temperature of the 
enzyme 44 ) in a final volume of 200 ul. The absorbance was recorded at 560 nm using a 
microtitre plate reader (Synergy HT, with Gen5 controlling software, Bioteck, USA). 
A control reaction was always carried out wherein 50 mM potassium phosphate 
buffer was added to the reaction mixture in lieu of the enzyme. Activity of SOD was 
expressed as units per mg protein. One unit of the enzyme activity was defined as the 
amount of enzyme required to inhibit nitro blue tetrazolium (NBT) reduction by 
50%. Protein concentration was determined by Bradford method 48 using bovine 
serum albumin as a standard reference. 

Assessing oligomerization status of the purified SODs. Oligomerization status was 
analyzed by electrophoresis under reducing and non-reducing conditions wherein 
equal quantities (1 jig) of the proteins were run on 15% sodium dodecyl sulfate- 
polyacrylamide (SDS-PAGE) 49 gel in the presence and absence of 2-mercaptoethanol 
(2-ME) 50 . Protein along with the sample loading buffer was boiled before loading onto 
SDS-PAGE. In a separate gel, the boiling step was omitted 50 and the protein separated 
on SDS-PAGE was stained for SOD activity 47 . 

Separation of monomer and dimer by size exclusion chromatography (SEC). 

Purified protein was filtered through 0.22 jj.m filter and subjected to SEC using a 
Superdex™75 10/300 GL(10 X 300-310 mm) (GE Healthcare, UK) connected to an 
Akta Prime™ plus (GE Healthcare, UK) system. The column was equilibrated with 
50 mM potassium phosphate buffer, pH 7.8, containing 40 mM KC1. Protein was 
applied in a volume of 500 ul. The column was operated at 20°C at a flow rate of 
0.5 ml min -1 . PrimeView 1M software was used for protein monitoring and 1.0 ml 
fractions were collected using an automated fraction collector. After SEC, fractions 
corresponding to monomeric and dimeric position were evaluated for protein 
estimation as described in earlier section. 

Effect of mutations on enzyme activity. Enzyme was autoclaved (12T C, 

1 . 1 kg/cm 2 , 20 min) and the activity of un-autoclaved and the autoclaved enzyme was 

determined as described in earlier section. Further, 300 ng of each un-autoclaved and 

autoclaved enzyme was run on a native PAGE using 10% separating and 

4% stacking gel 49 . After electrophoresis, gel was stained for activity 47 . 

Enzyme assay was performed at temperatures ranging between ( — ) 10 to ( + ) 80 'C. 
For assaying at sub-zero temperatures, 50% glycerol was added in the reaction mix- 
ture to avoid freezing and the assay was carried out by maintaining the temperature in 
a 96 well plate inside a deep freezer. SOD assay at temperatures ranging between 
( + )10 to ( + )80'"'C was performed in 0.5 ml tubes maintained at desired assay tem- 
peratures on a thermo-mixer (Eppendorf, Germany). Separate blank and control 
reactions were performed for all assay temperatures. 

Thermal inactivation assay. Thermal inactivation assay was performed' 8 with minor 
modifications. Purified enzyme was diluted to prepare 1-ml solution containing 
0.2 mg/ml protein in 50 mM potassium phosphate buffer, pH 7.8. For each 1 ml of 
the enzyme solution, a 100-ul of the initial sample was removed and stored on ice 
(0 min). Thereafter, the enzyme solution was heated in a water bath at 80°C. Every 
20 min, a 100 ul aliquot was removed till 160 min, stored on ice and assayed for 
residual SOD activity. The thermal inactivation data for WT and C95A was plotted 
for zero (residual activity versus time), first (natural logarithm of the residual activity 
versus time) and second order (reciprocal activity versus time) 18 . The rate constant fcj 
(min -1 ) for first-order thermal inactivation was determined 51 from the slope of the 
inactivation time course according to the "equation (1)", ln(A t /A 0 ) = — k^.t. Where 
A t is the residual activity that remains after heating the enzyme for time t, and A 0 is the 
initial enzyme activity before heating. The half-time of thermal inactivation (li/ 2 ) was 
determined 52 following the "equation (2)", t in = hi(2)/fcd- 



Limited proteolysis of enzymes. Enzymes were incubated with 1/20 w/w of trypsin, 
chymotrypsin or papain in 50 mM potassium phosphate buffer, pH 7.8, and 37 °C for 
defined time intervals. For chymotrypsin digestion, CaCl 2 was added to a final 
concentration of 20 mM in the buffer. For digestion with papain, cysteine was added 
to a final concentration of 5 mM in the buffer. Control samples were incubated under 
identical conditions but without the proteolytic enzyme. Aliquots were removed at 
defined time intervals and analyzed on 1 5% SDS-PAGE for protein integrity. Samples 
were also evaluated for activity on 10% native-PAGE. 

Circular dichroic spectroscopy and secondary structure estimation. Circular 
dichroic (CD) spectra were recorded on a JASCO J-815 spectropolarimeter equipped 
with a peltier thermostatic cell holder (model PTC- 423S/15; JASCO, Japan). 
Spectropolarimeter was pre -calibrated with 0.6% (w/v) aqueous ammonium 
( + )-camphor-10-sulphonate. The enzymes were dialyzed three times against 50 mM 
potassium phosphate buffer (pH 7.8). Far-UV (190-260 nm) CD spectra of enzyme 
(0.2 mg/ml) was recorded at 4 U C and 25°C in a quartz cuvette of 0.1 cm path length at 
a scan speed of 50 nm/min, with a 1 nm bandwidth, 1 nm data pitch and a Is 
response time. For thermal scans, the enzyme samples (0.2 mg/ml) were heated from 
0 to 100°C with a heating rate of 5°C/min controlled by a JASCO programmable 
peltier element. Far-UV CD spectrum was recorded every 10°C and the dichroic 
activity was monitored continuously. The results were expressed as mean residue 
ellipticity (0) by calculating mean residue weight per amino acid residue. Each 
spectrum was the mean of three accumulated scans. All values were corrected for 
solvent contributions. Molar ellipticity per residue (0), expressed in deg.cm 2 .dmol~\ 
was calculated using the "equation (3)", 0 MRE = MRW 0 o b s / 10d, where 0 is the 
ellipticity observed (mdeg), MRW is the mean residue weight per amino acid residue, 
c is protein concentration (mg/ml), and / is path length of the cuvette (cm). To 
estimate the proportions of secondary structures (a-helix, P-sheet, turns, random or 
unordered forms), the reference CD spectra obtained by Yang were used 53 . The 
software supplied by JASCO, Japan was used for analyzing the data. 

Molecular model-building. Models of three-dimensional structures of WT, C56A, 
C95A, and C145A were constructed using homology modeling approach. Structural 
based alignment was performed using Discovery Studio (DS) software package (M/S 
Accelrys, version 2.1, USA). X-ray crystallographic structure of WT (Protein Data 
Bank accession no. 2q2l) was used as a template structure. Root mean square 
deviation (RMSD), probability density function (PDF), total energy, molecular 
volume, and discrete optimized protein energy (DOPE) scores were calculated using 
the DS software. Model accuracy was determined by estimating the superimposition 
using "Verify 3D Profile" program of the DS. 
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